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^ ' Abstract 

O. 

CJ . A one dimensional disordered particle hopping rate asymmetric exclusion process (ASEP) with 

open boundaries and a random sequential dynamics is studied analytically. Combining the exact results 

^ ' of the steady states in the pure case with a perturbative mean field-like approach the broken particle-hole 

symmetry is highlighted and the phase diagram is studied in the parameter space {a, (3), where a and /3 



represent respectively the injection rate and the extraction rate of particles. The model displays, as in 



^+ . the pure case, high-density, low-density and maximum-current phases. All critical lines are determined 

QQ ' analytically showing that the high-density low-density first order phase transition occurs a,t a ^ f3. We 



show that the maximum-current phase extends its stability region as the disorder is increased and the 
usual 1/vt-decay of the density profile in this phase is universal. Assuming that some exact results for 
the disordered model on a ring hold for a system with open boundaries, we derive some analytical results 
for platoon phase transition within the low-density phase and we give an analytical expression of its 
corresponding critical injection rate a* . As it was observed numerically'- ^^\ we show that the quenched 
disorder induces a cusp in the current-density relation at maximum flow in a certain region of parameter 
space and determine the analytical expression of its slope. The results of numerical simulations we de- 
velop agree with the analytical ones. 
Key WordsiDisordered asymmetric exclusion model; steady state; boundary induced phase transitions; 

platoon phase. 
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1 Introduction 

It is well known that a set of equivalent problems including driven difFusion[l], Id turbulence [2], 
growth of interfaces [3] and directed polymers in a random medium in 1+1 dimension [4] is considered as 
the same problem described by the noisy Burgers equation [2] . All those systems exhibit a none equilibrium 
stationary states and show a wealth of interesting phenomena that don't occur in thermal equilibrium. 
The standard asymmetric exclusion process(ASEP) [1,5] is relevant to a variety of phenomena in physics[l] . 
It is the one which describes a driven lattice gaz with hard core exclusion. Under suitable scaling the 
evolution of its macroscopic density is governed by a noisy Burgers equation, which is equivalent to the 
KPZ equation[6], in Id. The ASEP with open boundary conditions is related to growth models with a 
defect or inhomogeneity[7] that introduces various types of phase transitions such as boundary induced 
phase transitions [8]. On the other hand it has a natural interpretation as a simplistic description of traffic 
flow on a single-lane high way and indeed forms the basis for more sophisticated traffic flow models [9]. 

The ASEP describes systems of interacting particles with hard core exclusion for double occupancy 
and with hopping rates differing for different directions. The bulk dynamics may be deterministic or 
stochastic. Open systems are coupled to reservoirs at both ends through stochastic boundaries. New 
methods for obtaining exact solutions for open Id ASEP systems have been derived for random-sequential 
update[10] and were solved exactly for any values of injecting and extracting rates at the boundaries [11]. 
At the same time a very elegant solution using a matrix product ansatz for the weights of the stationary 
states[12] was given. Therefore one can calculate bulk properties, phase diagrams, density profiles of 
boundary layers and interfaces between coexisting bulk phases as well as correlation functions. A domain 
wall (shock) approach was given[13] to understand quantitatively the system dynamics, stationary states, 
the nature and the location of phase transitions. The collective velocity and the domain wall velocity are 
the crucial ingredients of this picture. 

The different ways of updating sites are an essential part of the model. They affect the structure of 
phase diagram as well as the correlation function. The dynamics of updating may be applied in a 
random sequential order, parallel, i.e. fully synchronous for all sites or with sublattice parallel updating. 
For a detailed analytic results for ASEP we refer the reader to ref 14, 15 and references therein. The 
correlations are weakest for the random sequential updating, intermediate for ordered sequential and 
sublattice parallel updating and strongest for parallel updating. As a matter of fact, wealth of interesting 
exact results were obtained for random sequential and sublattice parallel dynamics, the fully parallel 
updating offers some difficulties. 

As far as a great deal of results are known about the ASEP with a deterministic jumps rates, a little is 
known analytically about systems with quenched disorder jump rates[16]. The random sequential version 
of the disordered model on a ring has been studied [17]. It was shown that a transition occurs between 



inhoinogcncous phase of low-density (jammed phase), where a traffic jam forms behind the slowest particle, 
and a high density congested phase (laminar phase), where all particles have to move more slowly than 
there preferred speed. It was shown that this transition persists in both dynamics parallel and ordered 
sequential updating and it is analogous to Bose condensation[18]. The disordered ASEP with open 
boundaries has been studied using numerical simulations[19]. A hopping parameter Ai was introduced 
to interpolate between fully parallel (At = 1) and random sequential (At -^ 0) dynamics. It was found 
that the phase diagram is similar to that obtained in the pure case a part the shift of the first order 
transition line describing the coexistence of the low and high-density phases. The numerical results in 
the maximum current phase are consistent with a power law-decay 1/vt of the density profile as in the 
pure case. 

In this paper we present an analytical study of a disordered ASEP with open boundary conditions. 
Particles jump in each time step to the right to vacant nearest neighbor sites with a probability p^ 
associated with particle /i . The probability of injected (extracted) particles at the right (left) boundary 
is a{P) . The model is defined in sec. 2 and studied analytically after developing the equation of the steady 
state probabilities[10] of a system of length N for a random sequential updating. Using a mean field-like 
approach we give an analytical solution for arbitrary a and f3. The disorder doesn't affect dramatically 
the topology of the phase diagram. We show in sec. 3 that the phase diagram exhibits three different 
regions: the low-density and high-density phases, which split into two regions characterized by different 
behavior of the density profile, and the maximum-current phase. The density profile in the various phases 
and the critical lines are obtained analytically and studied in sec. 4. Thus, an analytical expression of 
the shift induced by the broken particle-hole symmetry of the first order critical line between low and 
high density phases is obtained. The platoon phase transition is studied in sec. 5 where we discuss in 
detail the phase diagram and derive analytical expressions for bulk densities in low and high density 
phases. According to a power law distribution of particle hopping probabilities p^, we establish the 
analytical expression of the critical injection rate a* that corresponds to the platoon phase transition. 
We show that for the maximum current density, the current at this transition has gotten a negative slope 
whose analytical expression is calculated. To support our analytical results numerical simulations were 
performed in sec. 6. Our conclusion is given in sec 7. 

2 Model and Method 

We consider a one-dimensional lattice of N sites. Each site i, 1 < i < A^ , is either occupied by one 
particle or is empty. A configuration of the system is characterized by binary variables t^, where r^ = 
{Ti = 1) if site i is empty (full). A quenched random rate p^ is associated with each particle fi. The 
dynamics of the system is governed by the following rule: at each time step t ^ t + 1 one chooses at 



random an integer < i < N with a probability 1/{N +1), then the particle fi on site i hops to site 
(i + 1) if it is empty with a random probability p^ such that: 

Tt{t + 1) = 1 with probability p^Tj(f)Tj+i(t) + (1 -p^)Tj(t) 
Ti+i{t + 1) = 1 with probability Ti+i(t) +Pp(l - Ti+i{t))Ti{t) 

At the left boundary the site 1 remains occupied at time t + 1 if it was occupied at time i, and it gets 
occupied with probability a if it was empty at time t. Therefore: 

Ti(t + 1) = 1 with probability Ti(i) + a(l - Ti(i)) 
Ti(t + 1) = with probability (1 - q:)(1 - Ti(t)) 

Similarly, the site N remains empty at time i + 1 if it was empty at time t, and it gets empty with 
probability (3 if it was occupied at time t. So: 

TAr(i + 1) = 1 with probability (1 - l3)TN{t)) 

TAr(t + l) = with probability 1 - (1 - /3)(1 - TAr(t)) 

The steady state defined from the stochastic dynamics (l)-(3) is given in terms of the probabilities 
-PjvdTj}, {p^i}) of finding the specific configuration of particles represented by the set of the occupation 
numbers {ti} and jumping rate probabilities {p^} on the chain with N sites. Thus, following the dynamics 
rule described above it is easy to show that the probability PN{{Ti\, {p^j}) satisfies the relation: 

1 — a a 

PN{{T^}, {p^}) = J^^PN{{n}, {P^^) + ^^^Ti [Pn{0, {t,^i}, {p^}) + Pn(1, {n^i}, {pM 

+ ^^[^A^({n},K}) +Pa(t2 -ri)PAr(l,0,{n^i,2},{p^})] + ... + ^^ [/'wdrj, {p^}) 

1-/3 

+^^(1 - TN)[PN{{n^N}, 0, {p^}) + PN{{n^N}, 1, {pM^ (4) 

which can be written as: 

a(2ri - l)PAr(0,{T,^i},{p^}) + (t2 -Ti)pAPAr(l,0,{T,^i,2},{Pp}) + ... (5) 

+ {tn - TN-l)PvPN{{Ti^(N-l),N}, 1, 0, {pp}) + /3(1 - 2TAr)PAr({Ti^Ar}, 1, {p^}) = 0, 

where p\ and pi, are the hopping rate probabilities of particles A and v located respectively at site i — 1 
and site i ~ N — 1. 

It turns out to be more convenient to write the recursion relation with unnormalized probabilities 
fN{{n},{Pij}) related to /V({tj}, {p^}) by: 

PN{{n\, {p^,\) = ^-7T — r^ (6) 



where 

ZNi{p^})= E fNiinJAp,}) 

{n=0,l} 
In contrast to the pure case,i.e p^ = p, we are not able to find any recursion relation for fN{{Ti}, {Pfj.}) 
which allows to calculate exactly the average of any physical quantity as the occupation number, < r^ >. 
But, we will use a perturbative development based on a mean field- like approach to determine analytically 
such quantities. Indeed, if wc consider the unnormalized probability fN{{Ti},Pm) for the pure case, where 
Pm is the averaged value of p^ according to a common distribution (p{p) with a support on the the interval 
[c, 1], we may write: 



PN{{n},{p^}) 



^H% gVjv 



^w({Pm})' 

where Hn = ^"(/wdTj, {p^})), H^^ = ln(jN{{Ti},pyn)) and Vn = Hn - H%. Using a perturbative 
development of the potential Vjv that wc assume to be weak for jumping rate probabilities p^, we can 
determine the average occupation number < tj > as follow: 



< n >= 



Z^{Ti=0,l 



c^N e^N 



By introducing the partition function Zf^{pjn) for the pure case, the above equation can be written as: 



<n >- 



< r,e^" >o 



where the symbol < . >*' indicates that the average is made using the unnormalized probability /jv({Ti},pm 
Finally, using a perturbative development of the e^" we obtain 

< Ti >N^ (1- < Vn >")(< T, >° + < nVN >°). 

Knowing that for any real function f{xi, X2, ..., Xn) we may write, using the the differential operator dxi'- 

/(Ai, A2, ..., A„) == eEi=i ^*^"' f{xi,X2, ..., a:;„)/{j;,}=o, 

where Ai, A2, ..., A„ are real values, the potential VN(,{Ti}, {p^}) = '"•('fTT^rr^-^) is developed as follow: 

VN{{Ti], {p4) = eS" 1 -*S-- VV({xO, {Pm})/{..=0} (7) 

N 

= Y{e^^'^'VN{{xi\Ap,})Ux^=,} 

4=1 

Since r^ = 0, 1 the identity e'*''^' = 1 + Ti{e^ — 1) may be introduced in eq. 7 to obtain: 

N 

VAr({rJ,{p^}) =n[l + ^^(e^^* - l)]VN{{xi\,{P^.})l{..=^} (8) 






In order to evaluate the average < V/vdri}, {p^}) >'^ we use the dccouphng approximation < TiTjTk-.-Ti >'^~< 
Ti >°< Tj >°< Tk >° ... < T; >'' and neglect the higher order of the development, i.e (< Ti >° )™ <C 1 
for m > 1. Then, if we restrict ourselves to the first order we get from eq.8 

and then the average of the occupation number is given by: 



3 The density profile of the system 

In order to study the density profile of the system we compute from the recursion relation (6) the 
average occupation number < Ti >]\j such that: 



< Ti >Ar 



ZNilPt,}) 



(10) 



where 



TnAIP^^}) = X! TlfN{Tl,T2,...,TN,{Pf^}) (11) 

{Ti}=as 
The symbol < . > denotes the average over the configurations (ti,T2, ...,rjv) while the bar indicates the 
average on the quenched disorder. In order to analyze the density profile of the system, we will study its 
discrete derivative defined in ref 11. It is given for a set of values p^ by: 

TNJ+l{{p^,}) - TnAIPi^}) 



This quantity can be developed using eq. 10 and the pcrturbative development of eq. 9: 

N 



(12) 



t 



N 



{{Pf.})^ti{p^m)[l-VNiO,{pA) + J2<^^ >N (e^— l)VV({xJ,K})/{,,=0}] 



i=l 



where t^j^{pm) is the derivative of the density profile in the pure case for the averaged value p^ = p„ 
From the expression of the potential Vn{0, {Pfi}) we obtain 

fN{0,Pm) 



l-VW(0,{p^}) 



and 



(e^--l)VW({x,},{p4) 



.fN{0,{pA) 

fN{0,Pm.) d f fN{{xi},Pm) 



fNiO,{pA)dx, \fN{{Xi},{Pl^})J {x,= 



=0} 



Finally we get the expression of fj^HPf^})' 

4({P.}) = 4(p™)--^'^^''^"^ 



fNioApA) 



N 



I- 2^<Ti>j^ 



i=l 



dx, \fN{{xt}ApA)J {x. 



=0}. 



(13) 



In order to analyze the density profile of the system we will study the average value of its discrete 



derivative i^vdP/j}) obtained in eq. 13. However, to give the explicit form of t^rdP/^}) we have to 
determine the ones of < r^ >%, t^pf{pm) and the ratio fN{0,Pm)/ .fN{0, {P/j}) of unnormalized probabilities. 
The discrete derivative t^]^{pm) is obtained from eq. 12 since we may calculate the exact expression of 
ZN{Pm) and T^{pm) in the same manner as in ref 11. Thus, we obtain: 

4(p„0 = Z]^Hprn)il - — - -^){PpmyG^Z[j^_i{P/prn){aPrn)^-'G[iia/p^). 
Pm Pm 

It can be written more conveniently by using a simple rescaling a — > a/p^ and f3 -^ P/Pm as: 

4(Pm) = Z]^\pm)Fe{a/p^)FN-i{P/Pm) (14) 

with 

Fm{x) ^ x-^'-'G^^^ix) (15) 

and 

^^^"^ plJ^{l a/p^^ ~ /3/p„)(/?/p™)^+i(a/p™)^+i ^ ' 

ZN{Pm) can be given explicitly from the function Fj^{x): 

J^jv(/3/p^)--Fjv("/p^) for a^ 3 (v - 3) 

(a/p„)(l-a/p„)-(/3/p™)(l-/3/p„) ^T a ^ p , l,p^ p) ^^^^ 

~^Mf^ fora = /3^p„/2 

where F'{x) is the derivative with x. For the definition of the functions GfJ j^{x) we refer the reader to 
ref 11. 

On the other hand, the recursion relation of ,fN{{Ti},Pm) solution of eq. 5 can be constructed by the 
same way as in ref 10. We obtain: 



/Af(Tl,T2, ...,TN,Pm) = ap„iTNfN-liTl,T2, ...,TN,P,n) +/3Pm(l - TAr)...(l - Ti)/jV-l (ti , r2 , ...,TN,Pm) 

-al: 



J!/3 2_,(1 ~ TAr)...(l - Tfe+i)Tfe[/jV-l(Tl, ...,rfc-l,l,rfc+l, ...,TN,Pm) + fN-l{Tl, . .. , Tfc-l , 0, Tfc+l , ...,TN,Pm)], 



which leads to the unnormalized probability fN{0,Pm) = P{PPm)^^^- Although we are not able to give 
such recursion relation in the disordered case, we may derive from equation 4 the following recursion 
relations: 

/jv(0,{p^}) =/3p^/Ar-i(0, {p^}) and /jv(0, ..,0, 1, {p^}) = ap^/Ar-i(0, {p^}), 

where pi, denotes the hopping probability of the particle v located at site N —\. These recursion relations 
solve the steady state equation afN{0, {p^,}) — /3/jv(0, .., 0, 1, {p^}), obtained immediately from equation 
4, and easily lead to /Ar(0, {p^}) = PH^^^ l3p^. 



Up to the first order of the development in eq. 13, the discrete derivative of the average occupation is 

written as: 

N —, r- 

" " (18) 



Pr, 



By averaging over the disordered jumping rate probabihties it becomes: 



4(p{m}) = Z^^{p^)Ft{a/p„,)Fi_,{l3lp^) 



(19) 



where 



N 



F'n{^) 



vip) 



dp Fn{x) 



4 Density profile in the hydrodynamic limit and phase diagram 

The behavior of the density profile is discussed in the hydrodynamic limit iV ^ oo at large distances 



from both ends, i.e., £ >> 1 and r = (N — £) >> 1. The behavior of the average quantity t^j^{{p^}) allows 
us to localize the critical lines and then give the structure of the phase diagram. Using the asymptotic 
expression of FN{<j/pm) for large A^[ll]: 



l-2o-/p„ 



FN{<y/Pn 



^ if Cr < Prn/2 
y^(l-2cr/p„)2Ar3/2 if O- > Pm/2 



_ ) l(<T/p„)(l-a/p„)l« 
4" 



(20) 



the shape of the density profile is computed from eq(14-19). 

Since the method we use is a mean-field like approximation, it doesn't give the exact behavior near 
critical lines. But, it may be considered as a good approximation far from them as it is illustrated by 
the numerical results presented in what follows in sec. 6 (Fig 2). Thus, based on the approached density 
profile given in eq. 18 we will deduce the shape of the phase diagram from the divergence of scale lengths 
defined below. 

The high-density phase HDi is located within the region (] < a < Pmj'^- From the results obtained 
above the density profile decays exponentially such that: 



4({pp}) 



Pr, 



2a 



1- 



/^(Pm - ii) \ £/J 



where the length scale ^ is given by: 



r 



h 



a{pm - a] 
f3{Pm - /3) 



a{pm - a) 



fie) 



(21) 



(22) 



with/(c)=p„J^dp. 

c 

We remark that the length scale ^^^, which is rather different than the one obtained in the pure case, has 
gotten an additional factor that is responsible for the particle-hole symmetry breaking. The low-density 



phase LDi is bounded by the region defined by a < p„i/2 and (3 < Pm/'i with a < f3. As for HDi, the 



derivative of the density profile tj^{{p^}) may be obtained from eq. 14-19 such that: 



(23) 



where r — N — £ and the length scale (^ ^ — In 



a{pm~a) 



The current for high(low) density phase 



_/(c)/3(p„^/3)_ 

HDi(LDi) is given by j = Ppbuikij = a(l - Pbuik)) where pbuik = < tn >■ 

The coexistence line LDi ^^ HDi is determined by the divergence of the length scale £, defined in eq. 22. 

It is given by: 

^ - 2 U fie) J ^ ^ 

The high-density phase HD2 is located within the region defined by a > Pm/2 and f3 < Pm/'2- The 



behavior of the derivative i^ (p) changes to the form: 



^N 



V V'^iPrn - 2ay J 



where the length scale ^/3 is defined as: 



a' = In 



,P{Prn-0) 



Pn 



/(c) 



(25) 



(26) 



The current and the boundary values of the density are given by the same expressions as in the high- 
density phase HDi: < Tpf > — pbuik and < ti > = 1 — — pbuik- The profile of the low-density phase 
LD2(q: < Pm/'2 , /3 > Pm/'2) may be obtained using the same calculations as above. The density and the 
current are the same as in LDi phase but the density profile is given by: 






(27) 



where ^^1 — In 



4f{c 



a(pm—o:) 



The maximum-current phase is located in the region defined by a > Pm/'2 and (3 > Pm/'^- The 



derivative tj^{p) depends neither on a nor on /3. It takes the same expression as in the pure case 
independently of the disorder distribution: 



4({pm}) 



^ (-1 _ A^-3/2^-3/2 

4V^^ n' 



(28) 



However, the density approaches its bulk value pbuik — 1/2 as r ^" with the distance r = £ from the 
origin from above and r — N — £ from the boundary from below. The current depends on the average of 



the hopping rate parameter and it takes its maximal value j{{Pii}) = Pm/^- Therefore one obtains the 

boundary values < tm > — t^ and < ti > = 1 — ^. 

The phase transition between HD2 and the maximum-current phase MC, which is of second order, depends 



on the value of c for a given distribution ip{pfj). It occurs when the length scale ^^3 defined in eq. 26 
diverges and then is located at the critical value Pc'- 



P-n 



/3c = ^(l-v/l-l//W) (29) 

The density profile and the current should be continuous at the critical line. 

In order to determine the phase transition between LDi and LD2 we consider, for low values of a 
[a < ac), the transition between HD2 and MC phases on the line (3 ~ f3c- In the MC phase the bulk 
density pbuik ~ 1/2 and the way in which it is approached doesn't depend on (3, whereas in the high 
density phase HD2, /3 determines completely pbuik and the way the profile decays to it. This is due to the 
fact that for /3 > f3c, the particles close to the boundary don't block each other and are quickly extracted. 
As a result, any perturbation corresponding to a change in the extraction rate (3 doesn't spread into the 
system. As this description of the effect of /3 increasing beyond (3c on the transition from the high density 
phase IID2 to MC phase doesn't depend on the injection rate a at the origin, we conclude that the 
transition between low density LDi and LD2 phases should be located at the same critical line [3 — j3c- 
Based on the divergence of the correlation length ^q, (eq.27) one can find that the critical line LD2 ^^ MC 
is located at a = etc = /3c. It is of second order. The behavior of the profile at the critical line between 
high density phases HDi and HD2 may be determined using the same arguments as above based on the 
transition line LD2 ^^ MC. Consequently, the critical line between high density phases HDi and HD2 
should be at a = ac since the transition from the low density phase LD2 to MC phase doesn't depend on 
the extraction rate (3. Since all critical lines should meet at the critical point [ac, f3c) the coexistence line 
of LDi ^^ HDi described by the equation 24 should give the critical value f3c for a = otc- Unfortunately it 
is not the case. This discrepancy is due to the fact that our perturbative mean field-like approach doesn't 
give the exact value of the correlation length ^ as high correlations are neglected. As a consequence, the 
thermal average and the quenched disorder are, up to the first order of the development, decoupled (eq. 
18). We think that if we take into account of correlations we will obtain the exact value of ^ and then 
the discrepancy will be discarded. But, it is not obvious to develop any analytical result of the model 
using our approach without the decoupling approximation. In order to improve the expression of /3 that 
gives the correct critical line we remark that if we multiply the second term of eq. 24 by , ^ the 

critical value /3 = /3c is found for a ~ ac- We note that by defining a new correlation length: 



/3 {Pm - /3_2 „. X 



Q^'(Pm-"')' 

with the transformation B ~ 13 — ^(1 ^ ), a — , " and p^ = , ^"^ we recover 

^ ^ 2 V ^i+i//(c)^' Vi+V/(c) ^"^ V1+1//W 

that all lines meet at the critical point [ac, f3c). Such correlation length may be thought as resulting from 
the coupling of thermal average and quenched disorder since the new parameters a and (3 depend on 
the disorder through /(c). 

10 



On the other hand, (3 should vanish for a ^ while in the limit of the pure case, i.e /(c) -^ 1 we should 
obtain [3 ~ a. Consequently, a perturbative development to higher order should lead to the following 
expression of the coexistence critical line LDi ^^ HDi: 



1 / 2 t \\ 1/2 / 

Pm 1 Pm a[Pm-Ol)\ p,n a 




(3 = «(c) = ^ - , ^ - ^""; , ^' + £^ _ _ 1 1 _ (30) 

5 Discussion of results 

In this section we will discuss in detail the phase diagram, derive analytical expressions of bulk 
densities, pbuik, in high-density and low density phases and present some analytical results for platoon 
transition. 

5.1 Phase diagram and bulk densities 

The phase diagram (Fig.l) of the disordered ASEP presents the usual phases namely the high-density, 
low-density and maximum-current phases. The high-density phase is split into two regions. The argument 
of the logarithm of the length scale ^ in the first region (a < ac, [3 < (3c) depends on the distribution 







Pbuik = (1 )g(a,/3) 

Pm 




N 

ni 


(Pm\ 
\PJ 


N ^f / \ o 


f fN{x^,p^j)\ 
\fN{xuPm)) 



of the particle jumping rate through /(c) = p„i(l/p^) that is responsible for the particle-hole symmetry 
breaking. When S, diverges, a phase transition between high and low density phases occurs aX a ^ (3. 
The transition line which is of first order is given by eq. 48. Using the perturbative development of the 
bulk density pbuik — < tjv > we obtain, up to a second order, for high density region: 

(31) 

Pm 

where 

9{a,l3) 

^^, \ P„. / —r ^^ V Pu. I <JXi \iN\Xi,Vm.] I 

{2:.=0} 

Since we can not determine explicitly the expression of the function /Ar({a;i}, {p^}) we will give, based on 
some physical arguments, an approached form of the function g(a,/3) after having shown that it doesn't 
depend on the injection rate a. 

It was shown[17] that as the probability distribution of the gap (the number of vacant sites) in front of 
the /i*'' particle satisfies the stationary condition one provides p^a^ — constant — v, where a^ is the 
probability that the /i*'* particle has a vacant site in front of it and v is the average velocity of particles. 
Averaging over the quenched disorder we may obtain an approximated form of the velocity v{p) such that 
v{p) — pm'f^ — Pm(l — p)- Thus, the expression of the current of our disordered model may be written 
as j = pmp{l — p). The local density p{x,t) of the driven lattice gaz we study satisfies the continuity 
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equation g^' ' ''^ ' ' — 0. Consequently, the collective velocity may be given from the exact non 



equilibrium fluctuation dissipation tlieorem[13]: 



ycoll ^ 9j{p) 



dp 
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Figure 1: The fundamental diagram of the ASEP. The continuous lines present the analytical result of 
the critical lines while the squares show the numerical ones for L = 2000, c = 0.5 and n = 1. The dashed 
line(stars)localizes the platoon transition obtained analytically (numerically). The dotted lines delimit 
the different behaviors of low and high density regions. The inset represents an enlarging of the first 
critical line LDi-HDi. 



One finds X/^°" = p„(l — 2p). ft changes sign at p = f/2 where the current takes its maximal value 
3 ~ Pm/'^- fn high density regime and for a fixed j3, any small perturbation causes incoming particles to 
pile up behind the perturbation rather to spread into the bulk. As a result, the collective velocity of the 
center of mass of the perturbation is negative and thereafter the bulk density, pbuik , is independent on 
the injection rate a. Consequently, g{a,(3) = g{[3). On the other hand, one should recover the limiting 
cases namely that for c = 1 (pure case) g{j3) -^ 1 while in the limit /3 ^ the bulk density reaches its 
maximal value pbuik = 1- I3y increasing [3 for a fixed value of a > oic^ ptuik should decrease going to 
Pbuik = 1/2 at /? = /3c where a continuous phase transition to the MC phase occurs. Following these 
arguments the explicit expression of the bulk density in high density phase may be given by: 



(1 f^ 

Pbuik — t 

Pm 



f- 



I3{p,n - 2(5,) 

Wcipm - (3) 
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which may be simply written as: 

Pbuik = 1 - ^ (32) 

As for the pure case, the domain waU(shock)velocity is derived from the usual continuity equation: 

Vsh = (33) 

P+ - P- 

where jj-, p± are the bulk stationary state values of the current and the density in the left (-) and the 
right (+) parts of the domain wall. Since i± = PmP±(l ^ P±) we obtain 

Vsh == Pm(l - p+ - P-) (34) 

The expression of the domain wall velocity Vah changes its functional form depending on the phase 
transition lines where it should vanish. The bulk density in the low-density phase may be obtained by 
assuming that the domain wall picture[13] remains still valid for the disordered case. On the one hand, 
it should depend only on the injection rate a since any perturbation of the stationary state, for a fixed 
value of a, doesn't penetrate in the bulk. On the other hand, for the low-density/high-density domain 
wall (0|1) Vah should vanish at the first order critical line LDi ^^ HDi. As /9_|_ is given by eq. 32, we 
obtain the bulk density in the low-density phase pbuik — 2^' ^'^"'- ^^^ maximum current/high density 
domain wall (m|l) we have Vsh — ^ig — l)i(P- = 1/2)- It vanishes exactly at P — 13c- For the low 
density/maximum current domain wall (0|?7i) the domain wall velocity, Vgh — ^{^^^^ ^ 1)j vanishes at 
a ~ Oc, since K{a^ — a^. 

The shape of the phase diagram doesn't depend on the distribution of the particle jumping rate 
and MC phase gains more space by decreasing c since (Xc decreases. Indeed, it's easy to show that 
1 !i| ./(c) |< Pmlnc^^, which implies that /(c) > 1 and then ac -^ JITa for c << 1. 

It's worthwhile to note that the critical point (ac, (3c) where the three phases meet should be located in 
between e/2 and 1/2, which correspond respectively to the pure case p^ — c and p^ = 1. For the common 
distribution 

^(rt= (i!t)»+i (^'")" ^^^) 

we have Uc = (1 + n + c) n+2 ' which agrees with the value expected in ref 19. As it was argued 

numerically [19] the density in the MC phase decays algebraically and the deviation from its bulk value 
Pbuik — 1/2 decays asymptotically as Spi —\< Tg > ^pbuik |= a£^^'^. It is insensitive to the order of the 
maximum current. On the other hand, Hager et al[20] have conjunctured, using scaling arguments that 
this power law behavior holds for any one dimensional lattice gaz model. Based on the height difference 
correlation function of the one dimensional KPZ equation they have shown that the constant a, which 
was calculated for n > 1, is proportional to the compressibility. This result is independent on the shape 
of the distribution (p{p), the dynamics of updating[19,20], on the symmetry of the exclusion model[20] 
and neither on a nor /3. It is a universal behavior. 
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5.2 Platoon phase 

The exact current-density relation in the disordered ASEP can be obtained analytically in the hydro- 
dynamic limit N ^ (X) from an implicit equation for the stationary state particle velocity i;(p)[17,19]: 

P-(l + -J0-/p] (36) 

This relation holds for densities p* < p < 1 , where p* is the critical density for the onset of platoon 
formation[17,19]. For p < p* the overall speed will be set by the slowest particle and then v{p) = c while 
for p > p* , V is determined from eq. 36. Consequently, the value of p* is found by setting w = c in the 
rhs of this equation and the current-density relation becomes: 

M-h""' '"""' (37) 

[ pc li p< p* 

Since in the low-density phase we have Pbuik == K{a)/2ac the critical value a* for the platoon transition 
(The dashed line in fig. 1) is given by: 

1 

^^fj^^[p)dp (38) 

K[a*) J p-c 

c 

Consequently, the critical line separating the inhomogcneous phase (p < p*) from the congested phase 
(p > p*) is independent of (3 as it was shown numerically [19]. It depends only on the choice of the 
jumping rate probability distribution (p{p). According to the common distribution (35) we obtain 



(n + c) 



in agreement with the expression obtained in ref 19 since p* — '^„^". ' = "^ '^J . 

It was shown that the platoon phase transition is of second order, in the sense of ref 17, for < n < 1. 
In this regime the velocity v{p) is continuous at p ^ p* while it is discontinuous for n > 1. The current- 
density relation(eq. 37) leads to a shape with a quadratic maximum at pmax > P* i ^s it was shown in ref 
19. For n > \ the platoon transition becomes of first order and it is possible to give, within our approach, 
the value of c under which pmax — P* ■ As we have pmax — 1/2 we get for n > 1, c — jr^- 



From the expressions of the bulk density and the current given above, the derivative of j({p^}) with p is 
calculated at the platoon transition. We obtain for p* = Pmax, 

dj _J ^ ifp<P*-Pmax 

1 



-'max 



where i?(c) = p^ax J j^^^(pip)dp. 

c\ r\V-\i- o I T-k iV ( /-> 1 

2n(n-l)(l-c)2 
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According to the distribution 35 one obtain i){c) = 2n(n-i)(i-c)^ ■ ^^'^ easy to see that within the region for 



-'max- 



which p* = pmax one has l/'d{c) > c. Thus, the current presents a negative local slope for p > p* = p„ 
Consequently, the disorder induces a cusp at the maximum oi j{p) at the first order platoon transition. 
This result was observed in a previous numerical study [19]. 

6 Numerical simulations 

In order to support our analytical approach we have performed numerical simulations of the disordered 
ASEP on a lattice of sizes L = 2000 — 5 x 10'* and random initial distribution of particles. Once the 
system reaches the stationary state we compute the average < Q > of any physical quantity Q{t) over 
4 X 10^ — 10^ time steps. We note that, in contrast to the case of ring, a separate disorder average is not 
necessary in the case of open boundaries since new particles are constantly injected into the system. In 
what follows we will just give some numerical results that we compare with the analytical ones, namely 
the phase diagram, the density profile in the high density and MC phases and the platoon transition. 

The phase diagram(Fig. 1) is determined by computing the bulk density pbuik for different values of 
a and /3. For a fixed value of a < ac, where Oc is the critical value which separates the two high density 
regions HDi and HD2, ptuik decreases when increasing /3 until a critical value (3cr{a) where a first-order 
transition characterized by a discontinuity of phuik brings the system into the low-density phase where 
the density becomes independent of /3. The first-order transition high density-low density phase occurs 
at a critical line f3cr{(x) 7^ a. This is an effect of the disorder induced in the model which breaks the 
particle-hole symmetry. We note that the critical line obtained analytically is in good agreement with the 
numerical results. For a > ac the system exhibits, at a critical value /3c, a second-order phase transition 
to the MC phase where the bulk density and the current are independent of a and /3. The numerical 
and analytical values of /3c are in good agreement. For a fixed value of /3 > /3c the bulk density increases 
when increasing a until the critical value ac, where we find the MC phase. In order to check the validity 
of the analytical expressions of the density in the high density (eq. 32) and low density phases we have 
shown in fig. 2 the variation of pbuik in high(low) density with the extraction(injection)rate (3{a). We 
may consider that the analytical forms present a good fit of the numerical values especially for high values 
of c. 
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Figure 2: The variation of the bulk density vs /3 and a respectively in (a) high density and (b) low 
density regions for c = 0.75 and n — l.o presents the analytical expression and + denotes the numerical 
simulation result for L = 2000. 

The regions of the high density phase and low density phase exhibit different density profiles. In fig. 
3 we restrected oursclf to present the density profiles of HDi and HD2 phases. They present different 
behaviors far from the bulk from below as we have shown using our analytical approach. 
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Figure 3: The density profile of HDi (a = 0.2, /3 = 0.1) and HD2{a ^ 0.6, f3 = 0.1) phases. 

The density profile in the MC phase decays, within the mean field theory[8,21], as p{x) ^ x^/''"^^^ 
where m is an exponent that describes the behavior of J{p) — J{pmax) near the maximal density Pmax- 
The deterministic limit p^ ~ 1 oi the totally ASEP, within a fully parallel dynamics, for which m ~ I 
is the only exactly solved case and it is found[22] that the MC phase disappears altogether. In the dis- 
ordered ASEP the current-density relation depends on the choice of the distribution and the dynamics 
of updating. But, it seems that the density profile is not[19,20]. In fig. 4 we have presented the density 
profile in the MC phase for the random sequential updating. The data are consistent with a l/yZ-decay 
as we have shown analytically and they don't depend neither on n nor c. 
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Figure 4: A Log-Log plot of the density profile in the MC phase for L = 2000, n = 1, c = 0.5, a = 0.7 

and (3 — 0.9. The straight line indicates the l/vt-decay obtained analytically. 



In order to check the validity of our mean field-like theory and emphasize the analytical results 
obtained using this approach, we have shown in fig. 5 the variation of the critical value a* of the platoon 
transition according to n and c. It can be located by monitoring the variance of the head ways [17, 19] 



A"^ ^ <ul>f,-<Ua> 



Ai ^M 



(40) 



where < . >^ means the average over all the values of the gap u^ = .t^+i — .t^ — 1 (x^ is the position 
of the particle /i), as function of the injection rate a, i.e p, and system size. For a > a* (p > p*), A^ 
is independent of system size, while for a < a* (p < p*) it is dominated by the macroscopic gaps and 
acquires an L-dependence. 

In order to determine the value of a* , we compute the headway A^ for two system sizes namely L = 2000 
and L ~ 500. The data for the two system sizes coincide for a > a* but differ for a < a* . Since we 
restricted ourselves to high values of c, i.e c > 1/2, we remark, that the numerical data and the analytical 
expressions are in good agreement. The slight difference they may present for low values of c results from 
the fact that for such values the system presents high correlations which are neglected in our analytical 
approach. 
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Figure 5: The variation of a* vs (a) n (c = 0.5) and (b) c {n = l).o present the analytical result while + 
denotes the numerical data. The smooth curve represents the fit of the analytical results. 

7 Conclusion 

The effect of particle wise disorder on the phase diagram of the ASEP with open boundaries is studied 
analytically for a random sequential dynamics. The steady state has been established and its properties 
as the density profile the current and the bulk density were solved using the exact results of the pure case 
combined with a mean field-like approach. The shape of the phase diagram was found to be universal. 
It doesn't depend on the form of the distribution probability of particle jumping rate (p{p). The critical 
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line of the phase transition between the high-density phase and the low-density phase, which is of first 
order for any disorder distribution, was determined analytically. Because of the particle-hole symmetry 
breaking, we have shown that this critical line is located ai (3 — «(«) ^ a and its explicit form was given 
within our mean field-like approach. The phase transition between the high(low)-density phase and the 
maximum-current phase, which is of second order, is located at a critical value Pdoic ~ Pc)- We have 
shown, within our analytical approach that the value of ac depends only on the particle jumping rate 
distribution probability. Based on the fact that the transition between the maximal current phase and 
the high density phase(low density phase) doesn't depend on the injection(extraction)rate a(/3), we have 
deduced that the system exhibits a phase transition between high(low)-density phases HDi(LDi) and 
HD2(LD2) at the same critical value a — ac- The effect of the disorder doesn't change the nature of 
those critical lines. They remain of second order and they all meet at the critical point (ac, /3c — (Xc) that 
depends only on ip{p). As it was shown in a previous work[21], we have highlighted the universal behavior 
of the density profile in the maximum-current phase namely that the density of the system reaches its 
bulk value pbuik = 1/2 as l^^'"^ independently on the distribution (p{p) neither on a nor f3. 
Assuming that the analytical results obtained for the platoon phase transition occurring in the low-density 
phase for systems with periodic boundary conditions hold for systems with open boundaries[19], we have 
calculated the value of the injection rate a* at which the platoon phase transition occurs and we have 
shown that for a particular choice of the disorder distribution (p{p) (eq. 35) the transition is of first order 
for all values c < „ " . . Using the exact current-density relation (eq. 37) we have shown that the disorder 
induces a cusp at maximum flow in a certain region of parameter space as it was found numerically [19]. 
From our approach we were able to establish the analytical expression of the variation of the slope at 
maximal current. 

To check the validity of our analytical results obtained within our mean field-like theory, numerical 
simulations have been performed on a lattice of linear size. They show that the obtained results concerning 
the phase diagram and different quantities that are of interest are in good agreement with the analytical 
expressions. 

We think that some results obtained for a disordered model within a random sequential dynamics hold 
for the ordered sequential dynamics and the parallel dynamics. This point may be clarified in a future 
work. 
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